The Quantum— Mechanical Position Operator in Extended Systems 
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The position operator (defined within the Schrodinger representation in the standard way) becomes 
meaningless when periodic boundary conditions are adopted for the wavefunction, as usual in con- 
densed matter physics. We show how to define the position expectation value by means of a simple 
many-body operator acting on the wavefunction of the extended system. The relationships of the 
present findings to the Berry-phase theory of polarization are discussed. 



The position operator within the Schrodinger repre- 
sentation acts multiplying the wavefunction by the space 
coordinate. This is trivial, but only applies to the bound 
eigenstates of a finite system, which belong to the class 
of square-integrable wavefunctions. This is not the way 
how condensed matter theory works: almost invariably, 
one considers a large system within periodic boundary 
conditions (PBC), and the position operator (defined as 
usual) becomes then meaningless. For the sake of simplic- 
ity, most of this Letter will deal with the one-dimensional 
case. The Hilbert space of the single-particle wavefunc- 
tions is defined by the condition ip[x + L) — ip^x), where 
L is the imposed periodicity, chosen to be large with re- 
spect to atomic dimensions. An operator maps any vec- 
tor of the given space into another vector belonging to 
the same space: the multiplicative position operator x 
is not a legitimate operator when PBC are adopted for 
the state vectors, since x'ijj{x) is not a periodic function 
whenever is such. Of course, any periodic func- 

tion of x is a legitimate multiplicative operator: this is 
the case e.g. of the nuclear potential acting on the elec- 
trons. Since the position operator is ill defined, so is its 
expectation value, whose observable effects in condensed 
matter are related to macroscopic polarization. For the 
crystalline case, the long-standing problem of dielectric 
polarization has been solved a few years ago [|l-3||: po- 
larization is a manifestation of the Berry phase |4Q, i.e. 
is an observable which cannot be cast as the expectation 
value of any operator, being instead a gauge-invariant 
phase of the wavefunction. Here we find a different, and 
more fundamental, solution: we arrive indeed at defin- 
ing the expectation value of the position in an extended 
quantum system within PBC, where the operator enter- 
ing this definition is simple but rather peculiar. Amongst 
the most relevant features, the expectation value is de- 
fined modulo L, and the operator is no longer one-body: 
it acts as a genuine many-body operator on the peri- 
odic wavefunction of A'' electrons. In the crystalline case, 
the present result can be related to a discretized Berry 
phase, and sheds new light into the physical meaning of 
the latter. 

We study a system of A^ electrons in a segment of 
length L, and eventually the thermodynamic limit is 



taken: L oo, A oo, and N/L = uq constant. At 
any finite L the ground eigenfunction obeys PBC in each 
electronic variable separately: 



*o(a;i. 



,xn) ^ '^oixi, . . . ,x^ + L, . . . ,xn). (1) 



We assume the ground state nondegenerate, and we deal 
with insulating systems only: this means that the gap be- 
tween the ground eigenvalue and the excited ones remains 
finite for L — > oo. Since the spin variable is irrelevant to 
this problem, we omit it altogether and we deal with a 
system of spinless electrons. Our major goal is defining 
the expectation value of the electronic position (A), and 
to prove that our definition provides in the thermody- 
namic limit the physical macroscopic polarization of the 
sample. 

Before attacking the main problem, let us discuss the 
much simpler case where PBC are not chosen, and the 
A^-particle wavefunction (called <I>o in this case) goes to 
zero exponentially outside a bounded region of space. We 
may safely use the operator X = ^i; ^'^d define the 

position expectation value as usual: 



(A) - ($o|A|$o) - dxxn{x) 



(2) 



where n{x) is the one-particle density. The value of (A) 
scales with the system size, and the quantity of inter- 
est is indeed the dipole per unit length, which coincides 
with macroscopic polarization. The operator A, being 
the sum of identical operators acting on each electronic 
coordinate separately, is by definition a one-body opera- 
tor. The expectation value of this same operator cannot 
be evaluated if the wavefunction obeys PBC: in fact X 
does not commute with a translation by L, and therefore 
is not a legitimate operator in the Hilbert space defined 
by Eq. (0). 

We are now ready to state our main result, demon- 
strated in the following. When PBC are adopted, the 
position expectation value must be defined through: 

(A) - ^Imlog (^ole'^^l^o). (3) 

The expectation value (A) is thus defined only modulo L, 
hardly a surprising finding since the wavefunction is peri- 
odical. The operator occurring in Eq. (j|) is a legitimate 
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one in the Hilbert space of periodic wavefunctions; as a 
tradeoff it is no longer the sum of operators acting on each 
electronic coordinate separately, and is therefore a gen- 
uine many-body operator. At variance with Eq. (|^), the 
expectation value {X) within PBC cannot be expressed 
in terms of the one-particle density, not even in terms of 
the one-particle reduced density matrix. The A^-particle 
wavefunction is explicitly needed. 

What remains to be done is to prove that our defi- 
nition of (X) provides the relevant physical observable 
in the thermodynamic limit. Notice that L — > oo is a 
tricky limit, since the exponential operator in Eq. (|^) 
goes formally to the identity, but the size of the system 
and the number of electrons in the wavefunction increase 
with L Q . We will show that the electronic polarization 
(dipolc per unit length) is: 

Pel = hm ;^Im log (^-ole'^^l^o), (4) 

where e is the electron charge. It is expedient to intro- 
duce the family of Hamiltonians: 

1 ^ 

H{c^)^^Y.^P^~ha)^ + V, (5) 

i=l 

where a is a real constant, and V is the many-body po- 
tential. A parametric Hamiltonian of this kind was first 
introduced by W. Kohn many years ago 0|, and sub- 
sequently used by different authors [^,1). The ground 
eigenstate of H{0) is precisely |\E'o); more generally, the 
state vector e*""^|^o) fulfills the equation 

ff(a)e'"^|*o) =^oe'"^|*o), (6) 

with an a-independent Eq. This does not warrant that 
it is an eigenstate since PBC, Eq. (^, are not fulfilled, 
except in the special cases where a is a multiple of 2n/L. 

Since by hypothesis Eq is nondegenerate, e^'^^'^^^\"^o) 
is the ground eigenstate of H{2tt/L): we may then use 
perturbation theory to expand it to leading order in 1/L 
in terms of the eigenstates of H{0)- However, the 
standard formulas perform an arbitrary choice for the 
phase of the perturbed eigenstate; in the most general 
case we write instead: 

where P — X^iLi-Pi the momentum operator. It is 
important to realize that the perturbative expansion is 
a good approximation whenever L is much larger than a 
typical atomic dimension, while the number of electrons 
in the wavefunction and the system size are irrelevant. 

Replacement of Eq. @) into Eq. (||) shows that the 
phase 7i is a most fondamental quantity, since: 



iX) ^ ^, (8) 

but we have not yet related it to any observable. In 
order to prove that e{X) / L for large L is indeed the elec- 
tronic polarization — as anticipated in Eq. — it will be 
enough to show that its time derivative coincides with 
the adiabatic electrical current flowing through the sys- 
tem whenever the Hamiltonian contains a slowly varying 
time-dependent term. We start from 

dt"^ 2^ V(^^|e.^X|^^^ 

where |^o) is the time derivative of the instantaneous 
adiabatic eigenstate. Substituting now Eq. in Eq. (|^) 
the phase factor cancels out; to lowest order in 1/L we 
get 

e d , , iehsr^,- , , (^'JPIvI'o) , , 

Z^W-;;^E(*o|*.)^^^+c.c., (10) 

where c.c. indicates the complex conjugate. In Eq. (|l0| ) 
the j = term is omitted from the sum, since I^Po) can be 
taken as orthogonal to |^o) with no loss of generality; fur- 
themore we have exploited time-reversal symmetry, ow- 
ing to which all the adiabatic instantaneous eigenstates 
'i'j can be taken as real. 

This concludes our proof. In fact the right-hand mem- 
ber of Eq. ( p^ is the electronic current flowing through 
the system when the potential V is adiabatically varied, 
a well known expression due to Thouless: Eq. (2.5) in 
Ref. The rest of this Letter is devoted to an analysis 
of our major result, Eqs. (^) and (^), and of its relation- 
ship to previous work. 

The special case of iV = 1 corresponds to a lone quan- 
tum electron diluted in a large sample. The position ex- 
pectation value, Eq. can then be expressed in terms 
of the periodic density as: 

(X) = — Im log / dx e*T^?i(x). (11) 
27r Jo 

A similar expression has been previously used by a few 
authors in order to heuristically follow the adiabatic 
time evolution of a single quantum particle in a disor- 
dered condensed system within PBC. The case > 1 
is qualitatively different, in that — as stressed above — the 
operator used in Eq. (^ to define {X) is a genuine many- 
body one. This is particularly remarkable in view of the 
fact that the physical observable is an integrated current: 
the current is a typical one-body operator, as in fact is 
the operator P in the right-hand member of Eq. (|l^). 
The case of independent electrons is also worth comment- 
ing. In this special case the A^-particle wavefunction is 
uniquely determined by the one-body reduced density 
matrix p{x, x') (which is the projector over the set of the 
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occupied one-particle orbitals): therefore the expecta- 
tion value {X) is uniquely determined by p. 

For the crystalline case, macroscopic polarization is 
presently understood as a manifestation of the Berry 
phase [Q, both for independent electrons and for 
correlated electrons [^,1). The definition of Eq. (^) re- 
duces to the well established ones in the crystalline case: 
for finite L the present findings can be shown to be equiv- 
alent to a discretization of the line integral defining the 
Berry phase. In this Letter we provide an explicit proof 
for the independent-electrons case only: the correlated 
case is not much different. 

Suppose we have a crystalline system of lattice con- 
stant a, where we impose PBC over M linear cells: there 
are then M equally spaced Bloch vectors in the reciprocal 
cell [0,27r/a): 

= XT*' s = 0,l,...,M-L (12) 
Ma 

The size of the periodically repeated system is L = Ma. 
The one-body orbitals can be chosen to have the Bloch 
form: 

^qs;m{x + r) = e"'=^l/'9s,m(2:), (13) 

where r = Za is a lattice translation, and m is a band 
index. There are N/M occupied bands in the Slater de- 
terminant wavefunction, which we write as 

N/M M-1 

i*o)=An n (14) 

m— 1 s— 

where A is the antisymmetrizer. It is now expedient to 
define a new set of Bloch orbitals: 

i/'g.,m(a;) = e"'^^'0,,,™(x). (15) 

We then recast the expectation value of Eq. (^ , after a 
double change of sign, as: 

(X) = -^Imlog(*o|^'o), (16) 

where |'^o) is the Slater determinant of the ip^s. Accord- 
ing to a well known theorem, the overlap amongst two 
determinants is equal to the determinant of the overlap 
matrix amongst the orbitals: 

(X) ^ -^Im log det S, (17) 

ZTT 

where 

Ssm,s'm' = dx^pl^{x)e-'^^lljq^,.m'{x). (18) 

Jo 

Owing to the orthogonality properties of the Bloch func- 
tions, the overlap matrix elements vanish except when 



Qs' ^ Qs +2tt/L, that is s' = s-t-1. The NxN determi- 
nant can then be factorized into M small determinants: 

det 5= Yl detS{qs,qs+i), (19) 

s=0 

where — in order to make contact with previous literature 
— for the small overlap matrix we use the notation 

Sm,7n'iqs,qs+i) = / V'q^ „(x)e"'^^'0,^^i,™'(a;), 
Jo 

(20) 

and i/'gAf ,m(a;) = ipqo,rnix) is implicitly understood (so- 
called periodic gauge). Replacing Eq. ( po|) into Eq. (^) 
we get 

Pel = -7^ lim Ini log ]T det 5(g„ g.+i), (21) 

s=0 

which concludes our equivalence proof. Eq. ( pl| ) coincides 
in fact with the well known expression of the modern 
theory of polarization ||l|,^ , obtained by King-Smith and 
Vanderbilt by defining a (continuum) Berry phase as a 
line integral, and then discretizing it. 

The discretization of the Berry phase was originally 
introduced for purely computational purposes, and is 
in fact routinely used in first-principles calculations 
The alternate path followed here to arrive at the same 
result shows that the discretization has instead a very 
basic meaning of its own. Macroscopic polarization can 
be cast as the thermodynamic limit of an expression in- 
volving the expectation value of a relatively simple and 
physically meaningful many-body operator as in Eq. (|^). 
This operator "extracts" the Berry phase from the square 
modulus of the many-body wavefunction, which embeds 
the relevant information about the relative phases of the 
one-particle orbitals. 

Several previous findings about macroscopic polariza- 
tion |jl|-^ apply to the present formulation as well: we 
report here a few of them for the sake of completeness. 
One is not interested in defining an "absolute" polariza- 
tion: the measured bulk quantity is always the difference 
AP between two states of the given solid, connected by 
an adiabiatic transformation of the Hamiltonian: 

AP = APnuci + A/"ci = / dtJ{t), (22) 
Jo 

where J{t) is the total (nuclear -I- electronic) current flow- 
ing throught the sample while the potential V is adiabat- 
ically varied. Notice that in the adiabatic limit At goes 
to infinity and J(t) goes to zero. The quantity of inter- 
est AP can be evaluated as a two-point formula, using 
the initial and final states only: for the electronic term, 
one evaluates Eq. (^) with both the final and the ini- 
tial wavefunctions, and takes the difference. The result 
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is only defined modulo e; a similar indeterminacy ap- 
plies to the nuclear term APnuci- Nothing can be done 
about this ambiguity of the two-point formula, which ul- 
timately stems from Thouless's quantization of particle 
transport jl^ . There is of course no indeterminacy 
if one trades away the two-point formula and performs 
instead the time integral in Eq. (p^), using for the elec- 
tronic current the right-hand term of Eq. ([l^) . 

Generalization of Eqs. and (^) to the three- 
dimensional case requires some care. Since {X) is ex- 
tensive, within a naif approach the x-component of APd 
would be defined only modulo e/L^, which becomes van- 
ishingly small in the thermodynamic limit. Fortunately, 
the problem is less serious than this, and the drawback 
is easily eliminated by adapting to the present formu- 
lation a major finding from Ref. j^. Suppose the sys- 
tem is crystalline with a simple cubic lattice of constant 
a. Then — upon exploiting the lattice periodicity of the 
one-particle density — it can be shown that the two-point 
formula provides each component of APd with an inde- 
terminacy of e/a^, which is no serious drawback. Such 
indeterminacy has nothing to do with electron correla- 
tion, and not even with quantum mechanics: a similar 
indeterminacy is also present in the classical nuclear term 
APnuci whenever this term is evaluated as a two-point 
formula . If the system is noncrystalline, then a large 
"supercell" is needed to reproduce the disorder, and only 
small polarization differences are accessible via the two- 
point formula. Again, this looks like a fundamental con- 
sequence of Thouless's quantization of particle transport 

Em- 

The modern viewpoint about macroscopic polarization 
has even spawned a critical rethinking of density- 
functional theory in extended systems. The debate 
started in 1995 with a paper by Gonze, Ghosez, and 
Godby jl^, and continues these days [|4). Although the 
subject is clearly outside the scope of this Letter, I point 
out that the present main achievement — namely, defin- 
ing the polarization of a many-electron system by means 
of an expectation value — could possibly help in further 
clarifying the matter. 
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